Fractional Quantum Hall Effect in Optical Lattices 
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We analyze a recently proposed method to create fractional quantum Hall (FQH) states of atoms 
confined in optical lattices [A. S0rensen et al, Phys. Rev. Lett. 94 086803 (2005)]. Extending the 
previous work, we investigate conditions under which the FQH effect can be achieved for bosons on 
a lattice with an effective magnetic field and finite onsite interaction. Furthermore, we characterize 
the ground state in such systems by calculating Chern numbers which can provide direct signatures 
of topological order and explore regimes where the characterization in terms of wavefunction overlap 
fails. We also discuss various issues which are relevant for the practical realization of such FQH 
states with ultra cold atoms in an optical lattice, including the presence of the long-range dipole 
interaction which can improve the energy gap and stabilize the ground state. We also investigate 
a new detection technique based on Bragg spectroscopy to probe these system in an experimental 
realization. 
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INTRODUCTION 



With recent advances in the field of ultra-cold atomic 
gases, trapped Bose-Einstein condensates (EEC's) have 
become an important system to study many-body physics 
such as quantum phase transitions. In particular, the 
ability to dynamically control the lattice structure and 
the strength of interaction as well as the absence of im- 
purities in EEC's confined in optical lattices, have led to 
the recent observation of the superfluid to Mott-insulator 
transition [U, i, i i, At the same time, there has 
been a tremendous interest in studying rotating EEC's 
in harmonic traps; at sufficient rotation an Abrikosov lat- 
tice of quantized vortices has been observed ,6] and the 
realization of strongly correlated quantum states simi- 
lar to the fractional quantum Hall states has been pre- 
dicted to occur at higher rotation rates [Z, 9] ■ In these 
proposals, the rotation can play the role of an effective 
magnetic field for the neutral atoms and in analogy with 
electrons, the atoms may enter into a state described 
by the Laughlin wavefunction, which was introduced to 
describe the fractional quantum Hall effect. While this 
approach yields a stable ground state separated from all 
excited states by an energy gap, in practice this gap is 
rather small because of the weak interactions among the 
particles in the magnetic traps typically used. In optical 
lattices, the interaction energies are much larger because 
the atoms are confined in a much smaller volume, and 
the realization of the fractional quantum Hall effect in 
optical lattices could therefore lead to a much higher en- 
ergy gap and be much more robust. In a recent paper 
, it was shown that it is indeed possible to realize the 
fractional quantum Hall effect in an optical lattice, and 
that the energy gap achieved in this situation is a frac- 



* Electronic address: ' hafezi@fas. harvard. edul 



tion of the tunneling energy, which can be considerably 
larger than the typical energy scales in a magnetic trap. 

In addition to being an interesting system in its own 
right, the fractional quantum Hall effect, is also inter- 
esting from the point of view of topological quantum 
computation [ll[ . In these schemes quantum states with 
fractional statistics can potentially perform fault tolerant 
quantum computation. So far, there has been no direct 
experimental observation of fractional statistics although 
some signatures has been observed in electron interfer- 
ometer experiments (T2| . [isj . Strongly correlated quan- 
tum gases can be a good alternative where the systems 
are more controllable and impurities are not present. 
Therefore, realization of fractional quantum Hall states 
in atomic gases can be a promising resource for topolog- 
ical quantum computation in the future. 

As noted above the FQH effect can be realized by sim- 
ply rotating and cooling atoms confined in a harmonic 
trap. In this situation, it can be shown that the Laugh- 
lin wavefunction exactly describes the ground state of 
the many body system In an optical lattices, on 

the other hand, there are a number of questions which 
need to be addressed. First of all, it is unclear to which 
extent the lattices modifies the fractional quantum Hall 
physics. For a single particle, the lattice modifies the 
energy levels from being simple Landau levels into the 
fractal structure known as the Hofstadter butterfly fTB\. 
In the regime where the magnetic flux going through each 
lattice a is small, one expects that this will be of minor 
importance and in Ref. [lO| it was argued that the frac- 
tional quantum Hall physics persists until a < 0.3. In 
this paper, we extent and quantify predictions carried 
out in Ref. [l^ . Whereas Ref. only considered the 
effect of infinite onsite interaction, we extent the anal- 
ysis to finite interactions. Furthermore, where Ref. [13 
mainly argued that the ground state of the atoms was in 
a quantum Hall state by considering the overlap of the 
ground state found by numerical diagonalization with the 
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Laughlin wavefunction, we provide further evidence for 
this claim by characterizing the topological order of the 
system by calculating Chern numbers. These calcula- 
tions thus characterize the order in the system even for 
parameter regimes where the overlap with the Laughlin 
wavefunction is decreased by the lattice structure. 

In addition to considering these fundamental features 
of the FQH states on a lattice, which are applicable re- 
gardless of the system being used to realize the effect, we 
also study a number of questions which are of particu- 
lar interest to experimental efforts towards realizing the 
effect with atoms in optical lattices. In particular, we 
show that adding dipole interactions between the atoms 
can be used to increase the energy gap and stabilize the 
ground state. Furthermore, we study Bragg spectroscopy 
of atoms in the lattice and show that this is a viable 
method to identify the quantum Hall states created in 
an experiment, and we discuss a new method to generate 
an effective magnetic field for the neutral atoms in the 
lattice. 

The paper is organized as follows: In Sec. [2l we study 
the system with finite onsite interaction. In Sec. [3] we 
introduce Chern numbers to characterize the topological 
order of the system. The effect of the dipole-dipole inter- 
action has been elaborated in Sec. 14.11 Sec. 14.21 studies 
the case of = 1/4 fiUing factor. Sec. [S]is dedicated to 
explore Bragg spectroscopy of the system. Sec.[S]outlines 
a new approach for generating the type of the Hamilto- 
nian studied in this paper. 



QUANTUM HALL STATE OF BOSONS ON A 
LATTICE 



2.1. The Model 

The fractional quantum Hall effect occurs for electrons 
confined in a two dimensional plane under the presences 
of a perpendicular strong magnetic field. If N is the 
number of electrons in the system and is the num- 
ber of magnetic fluxes measured in units of the quantum 
magnetic flux $o = h/e, then depending on the filling 
factor 1/ = N/N^ the ground state of the system can 
form highly entangled states exhibiting a rich behaviors, 
such as incompressibility, charge density waves, and any- 
onic excitations with fractional statistics. In particular, 
when 1/ — 1/to, where m is an integer, the ground state 
of the system is an incompressible quantum liquid which 
is protected by an energy gap from all other states, and 
in the Landau gauge is well described by the Laughlin 
wavefunction I16l: 



N N 

^{z,,z2, = Y[{z, ~ z^r n (1) 

j>k j=l 

where the integer m should be odd in order to meet the 
antisymmetrization requirement for fermions. 



Although the fractional quantum Hall effect occurs for 
fermions (electrons), bosonic systems with repulsive in- 
teractions can exhibit similar behaviors. In particular, 
the Laughlin states with even m correspond to bosons. 
In this article, we study bosons since the experimen- 
tal implementation are more advanced for the ultra-cold 
bosonic systems. We study a system of atoms confined in 
a 2D lattice which can be described by the Bose-Hubbard 
model [17] with the Peierls substitution jTsl. [Tsj. 
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where J is the hopping energy between two neighboring 
sites, U is the onsite interaction energy, and 27ra is the 
phase acquired by a particle going around a plaquette. 
This Hamiltonian is equivalent to the Hamiltonian of a 
U(l) gauge field (transverse magnetic field) on a square 
lattice. More precisely, the non-interacting part can be 
written as 



- J ^ a|aj exp 

<ij> 



2m 

$0 



A-dl 



(3) 



where A is the vector potential for a uniform magnetic 
field and the path of the integral is chosen to be a straight 
line between two neighboring sites. In the symmetric 
gauge, the vector potential is written as A = -f (— y, x, 0). 
Hence, a will be the amount of magnetic flux going 
through one plaquette. 

While the Hamiltonian in Eq. ^ occurs naturally for 
charged particles in a magnetic field, the realization of a 
similar Hamiltonian for neutral particles is not straight- 
forward. As we discuss in Sec. [6] this may be achieved 
in a rotating harmonic trap, and this have been very 
successfully used in a number of experiments in mag- 
netic traps [l9l . [20| , but the situation is more complicated 
for an optical lattice. However, there has been a num- 
ber of proposals for lattice realization of a magnetic field 
10, 21, 22], and recently it has been realized experimen- 
tally 23] . Popp et al. 2M have studied the realization of 
fractional Hall states for a few particles in individual lat- 
tice sites. A new approach for rotating the entire optical 
lattice is discussed in Sec. IH] The essence of the above 
Hamiltonian is a non-zero phase that a particle should 
acquire when it goes around a plaquette. This phase 
can be obtained for example by alternating the tunnel- 
ing and adding a quadratic superlattice potential |,10j or 
by simply rotating the lattice (Sec [6]). The advantage of 
confining ultra-cold gases in an optical lattice is to en- 
hance the interaction between atoms which consequently 
result in a higher energy gap comparing to harmonic trap 
proposals (e.g. Ref. |7|). This enhancement in the energy 
gap of the excitation spectrum can alleviate some of the 
challenges for experimental realization of the quantum 
Hall state for ultra-cold atoms. 
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2.2. Energy spectrum and overlap calculations 

In order to approximate a large system, we study the 
system with periodic boundary condition, i.e. on a torus, 
where the topological properties of the system is best 
manifested. 

There are two energy scales for the system: the first 
is the magnetic tunneling term, Ja, which is related to 
the cyclotron energy in the continuum limit Hujc — 47r Ja 
and the second is the onsite interaction energy U. Exper- 
imentally, the ratio between these two energy scales can 
be varied by varying the latticepotential height [l|, [13 
or by Feshbach resonances [1, [2^, [26i|. Let us first as- 
sume that we are in the continuum limit where a <C 1, 
i.e. the flux through each plaquette is a small fraction of 
a flux quantum. A determining factor for describing the 
system is the filling factor i' — N/N^, and in this study 
we mainly focus on the case of v = 1/2, since this will be 
the most experimentally accessible regime. 

We restrict ourself to the simplest boundary conditions 
for the single particle states ts{L)ip{xs,ys) = i'ixsTUs), 
where ts{L) is the magnetic translation operator which 
translates the single particle states 'ip{xs,ys) around the 
torus. The definition and detailed discussion of the 
boundary conditions will be elaborated in Section [31 The 
discussed quantities in this section, such as energy spec- 
trum, gap and overlap, do not depend on the boundary 
condition angles (this is also verified by our numerical 
calculation). 

In the continuum case, for the filling fraction v = 1/2, 
the Laughlin state in Landau gauge {A = —Byx) is given 
by Eq. ^ with m = 2. The generalization of the Laugh- 
lin wavefunction to a torus takes the form 1271 



'^{zi,Z2, Zpf) = frel{zi, ^2, ZN)Fcm{Z)e 



(4) 



where frei is the relative part of the wave function and 
is invariant under collective shifts of all z^'s by the same 
amount, and Fcm{Z) is related to the motion of the center 
of mass and is only a function of Z = Tli'^i- 1^°^ ^ 
system on a torus of the size {L^ x Ly), we write the 
wavefunction with the help of theta functions, which are 
the proper oscillatory periodic functions and are defined 
as d[l]{z\T) = ^^e*'^^("+'')'+2'^'("+'')(^+^) where the 
sum is over all integers. For the relative part we have. 



i<j 



(5) 



According to a general symmetry argument by Hal- 
dane 28j . the center of mass wave function Fcm{Z) is 
two- fold degenerate for the case of = 1/2, and is given 

by 



1/2 + (N^ - 2)/4 



^-^\2^^] (6) 



where Z = 0, 1 refers to the two degenerate ground states. 
This degeneracy in the continuum limit is due to the 
translational symmetry of the ground state on the torus, 
and the same argument can be applied to a lattice system 
when the magnetic length is much larger than the lattice 
spacing a <C 1. For higher magnetic filed, the lattice 
structure becomes more pronounced. However, in our 
numerical calculation for a moderate magnetic field a < 
0.4, we observe a two-fold degeneracy ground state well 
separated from the excited state by an energy gap. We 
return to the discussion of the ground state degeneracy 
in Sec. [31 



In the continuum limit a <C 1, the Laughlin wavefunc- 
tion is the exact ground state of the many body system 
with a short range interaction 0, [U, . The reason is 
that the ground state is composed entirely of states in the 
lowest Landau level which minimizes the magnetic part 
of the Hamiltonian, the first term in Eq. 0. The expec- 
tation value of the interacting part of the Hamiltonian, 
i.e. the second term in Eq. ([2]), for the Laughlin state is 
zero regardless of the strength of the interaction, since it 
vanishes when the particles are at the same position. 



To study the system with a non-vanishing a, we 
have performed a direct numerical diagonalization of the 
Hamiltonian for a small number of particles. Since we are 
dealing with identical particles, the states in the Hilbert 
space can be labeled by specifying the number of parti- 
cles at each of the lattice sites. In the hard-core limit, 
only one particle is permitted on each lattice site, there- 
fore for N particles on a lattice with the number of sites 
equal to {Nx = Lx/a,Ny = Ly/a)^ where a is the unit 
lattice side, the Hilbert space size is given by the com- 



bination 



On the other hand, in 



' ; Af!(Af,JV„-A')!- 

case of finite onsite interaction, the particles can be on 
top of each other, so the Hilbert space is bigger and is 
given by the combination ( " + " ^ ) . In our simu- 
lations the dimension of the Hilbert can be raised up 
to ~ 4 • 10^ and the Hamiltonian is constructed in the 
configuration space by taking into account the tunneling 
and interacting terms. The tunneling term is written in 
the symmetric gauge, and we make sure that the phase 
acquired around a plaquette is equal to 27ra, and that 
the generalized magnetic boundary condition is satisfied 
when the particles tunnels over the edge of the lattice 
[to be discussed in Sec. [31 c.f. Eq. (fTTj) ]. By diagonaliz- 
ing the Hamilton, we find the two-fold degenerate ground 
state energy which is separated by an energy gap from 
the excited states and the corresponding wavefunction in 
the configuration space. The Lauhglin wave function Q 
can also be written in the configuration space by simply 
evaluating the Laughlin wave function at discrete points, 
and therefore we can compared the overlap of these two 
dimensional subspaces. 
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FIG. 1: (color online) (a) The overlap of the ground state 
with the Laughlin wavefunction. For small a the Laughlin 
wavefunction is a good description of the ground state for 
positive interaction strengths. The inset shows the same re- 
sult of small U. (b) The energy gap for N/N^ = 1/2 as a 
function of interaction U/J from attractive to repulsive. For 
a fixed a, the behavior does not depend on the number of 
atoms. The inset define the particle numbers, lattice sizes, 
and symbols for both parts (a) and (b). 
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FIG. 2: (color online) (a) The energy gap as a function of all 
and a for a fixed number of atoms (N=4). The gap is cal- 
culated for the parameters marked with dots and the surface 
is an extrapolation between the points, (b) Linear scaling of 
the energy gap with all for U <^ J,C( < 0.2. The results are 
shown for N = 2(0), N = 3(*), iV = 4(-f) and iV = 5(V). 
The gap disappears for non-interaction system, and increases 
with increasing interaction strength (oc all) and eventually 
saturate to the value in the hardcore limit. 



2.3. Results with the finite onsite interaction 

The energy gap above the ground state and the ground 
state overlap with the Laughlin wavefunction for the case 
of 1/ = 1/2 in a dilute lattice a < 0.2, are depicted in 
figure. [TJ The Laughlin wavefunction remains a good 
description of the ground state even if the strength of 
the repulsive interaction tends to zero (Fig.[T]a). Below, 
we discuss different limits: 

First, we consider U > , U ^ Ja: If the inter- 
action energy scale U is much larger than the magnetic 
one (Ja), all low energy states lie in a manifold, where 
the highest occupation number for each site is one, i.e. 
this corresponds to the hard-core limit. The ground state 
is the Laughlin state and the excited states are various 
mixtures of Landau states. The ground state is two-fold 



degenerate ^2§\ and the gap reaches the value in for the 
hard-core limit at large U > 3J/a, as shown in Fig. [2] 
(a). In this limit the gap only depends on the tunneling 
J and flux a, and the gap is a fraction of J. These results 
are consistent with the previous work in Ref. [To| . 

Secondly, we consider |C/| ^ Ja. In this regime, the 
magnetic energy scale {J a) is much larger than the inter- 
action energy scale U. For repulsive regime (J7 > 0), the 
ground state is the Laughlin state and the gap increases 
linearly with all, as shown in Fig.[2]b. 

Thirdly, we study U — where the interaction is ab- 
sent and the ground state becomes highly degenerate. 
For a single particle on a lattice, the spectrum is the fa- 
mous Hostadfer's butterfly [l^, while in the continuum 
limit a ^ 1, the ground state is the lowest Landau level 
(LLL). The single particle degeneracy of the LLL is the 
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number of fluxes going through the surface, N^. So in 
the case of N bosons, the lowest energy is obtained by 
putting N bosons in levels. Therefore, the many- 
body ground state's degeneracy should be: ^ " jj, ^ ' 
For example, 3 bosons and 6 fluxes gives a 336-fold de- 
generacy in the non-interacting ground state. 

If we increase the amount of phase (flux) per plaque- 
tte (a) we are no longer in the continuum limit. The 
Landau level degeneracy will be replace by ^^^^ where 
a = r/ s is the amount of flux per plaquette and r and s 
are coprime [30| . Then, the many-body degeneracy will 

be: ( " i~ ' ) ■ 

Fourthly, we consider J7 < , J7 ^ Ja: when U is 
negative (i.e. attractive interaction) in the limit strong 
interaction regime, the ground state of the system will 
become a pillar state. In a pillar state, all bosons in 
the system condensate into a single site. Therefore, the 
degeneracy of the ground state is x Ny and the ground 
state manifold can be spanned as. 
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K) 



vac) 



(7) 



These states will very fragile and susceptible to collapse 
p. 

In a lattice, it is also possible to realize the fractional 
quantum Hall states for attractive interaction in the limit 
when |J7| 3> J a. Assume that the occupation number of 
each site is either zero or one. Since the attraction en- 
ergy U is very high and there is no channel into which 
the system can dissipate its energy, the probability for a 
boson to hop to a site where there is already a boson is in- 
finitcsimally small. Therefore, the high energy attraction 
will induce an effective hard-core constraint in the case 
of ultra-cold system. The energy of these state should 
be exactly equal to their hard-core ground state counter- 
parts, since the interaction expectation value of the inter- 
action energy is zero for the Laughlin state. The numeri- 
cal simulation shows that these two degenerate states in- 
deed have a good overlap with the Laughlin wavefunction 
similar to their repulsive hard-core counterparts and also 
their energies are equal to the hard-core ground state. 
These states are very similar to repulsively bound atom 
pairs in an optical lattice which have recently been ex- 
perimentally observed [3] . 

So far we have mainly considered a dilute lattice a < 
0.2, where the difference between a lattice and the con- 
tinuum is very limited. We shall now begin to inves- 
tigate what happens for larger values of a, where the 
effect of the lattice plays a significant role. Fig. [3l shows 
the ground state overlap with the Laughlin wavefunc- 
tion as a function of the strength of magnetic flux a and 
the strength of the onsite interaction U . As a increases 
the Laughlin overlap is no longer an exact description 
of the system since the lattice behavior of the system is 
more pronounced comparing to the continuum case. This 
behavior doesn't depend significantly on the number of 



particles for the limited number of particles that we have 
investigated iV < 5. We have, however, not made any 
modification to the Laughlin wave function to take into 
account the underlying lattice, and from the calculations 
presented so far, it is unclear whether the decreasing wave 
function overlap represents a change in the nature of the 
ground state, or whether it is just caused by a modifica- 
tion to the Laughlin wave function due the difference be- 
tween the continuum and the lattice. To investigate this, 
in the next Section, we provide a better characterization 
of the ground state in terms of Chern numbers, which 
shows that the same topological order is still present in 
the system for higher values of a. 

As a summary, we observe that the Laughlin wavefunc- 
tion is a good description for the case of dilute lattice 
{a ^ 1), regardless of the strength of the onsite inter- 
action. However, the protective gap of the ground state 





FIG. 3: (color online) Ground state overlap with the Laughlin 
wavefunction. (a) and (b) are for 3 and 4 atoms on a lattice, 
respectively, a is varied by changing the size of the lattice 
(the size in the two orthogonal directions differ at most by 
unity). The Laughlin state ceases to be a good description of 
the system as the lattice nature of the system becomes more 
apparent a > 0.25. The overlap is only calculated at the posi- 
tions shown with dots and the color coding is an extrapolation 
between the points. 
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becomes smaller for weaker values of interaction and in 
the perturbative regime U <^ J is proportional to all for 
a < 0.2. 



3. CHERN NUMBER AND TOPOLOGICAL 
INVARIANCE 

3.1. Chern number as a probe of topological order 

In the theory of quantum Hall effect, it is well un- 
derstood that the conductance quantization, is due to 
the existence of certain topological invariants, so called 
Chern numbers. The topological invariance was first in- 
troduced by Avron et al. 32] in the context of Thouless 
et al. (TKNdN)'s original theory [s^ about quantiza- 
tion of the conductance. TKNdN in their seminal work, 
showed that the Hall conductance calculated from the 
Kubo formula can be expressed into an integral over the 
magnetic Brillouin zone, which shows the quantization 
explicitly. The original paper of TKNdN deals with the 
single-particle problem and Bloch waves which can not be 
generalized to topological invariance. The generalization 
to many-body systems has been done by Niu et al.^^ 
and also Tao et a/.js^, by manipulating the phases de- 
scribing the closed boundary conditions on a torus (i.e. 
twist angles) , both for the integer and the fractional Hall 
systems. These twist angles come from natural general- 
ization of the closed boundary condition. 

To clarify the origin of these phases, we start with 
a single particle picture. A single particle with charge 
{q) on a torus of the size {Lx,Ly) in the presence of a 
magnetic field B perpendicular to the torus surface, is 
described by the Hamiltonian 



1 

2m 



-ih 



d_ 

dx 



d_ 

"dy 



where A is the corresponding vector potential (^f — 
= B). This Hamiltonian is invariant under the mag- 
netic translation. 



(8) 



is (a) = e 



(9) 



where a is a vector in the plane, and k'" is the pseudo- 
momentum, defined by 



d 

K = - lAx - qBy 

d 

ky — ~ lAy + qBx 



(10) 



The generalized boundary condition on a torus is given 
by the single-particle translation 



ts{Lxx)'ip{xs,ys) 
ts{Lyy)ip{xs,ys) 



''4>{xs,ys) 
'ip{xs,ys) 



where 0i and 62 are twist angles of the boundary. The ori- 
gin of these phases can be understood by noting that the 
periodic boundary conditions corresponds to the torus in 
Fig. m^a). The magnetic flux through the surface arises 
from the field perpendicular to the surface of the torus. 
However, in addition to this fiux, there may also be fluxes 
due to a magnetic field residing inside the torus or pass- 
ing through the torus hole, and it is these extra fluxes 
which give rise to the phases. The extra free angles are 
all the same for all particles and all states in the Hilbert 
space, and their time derivative can be related to the 
voltage drops across the Hall device in two dimensions. 

The eigenstates of the Hamiltonian, including the 
ground state will be a function of these boundary angles 
(01, 62). By defining some integral form of this eigen- 
state, one can introduce quantities, that do not depend 
on the details of the system, but reveal general topolog- 
ical features of the eigenstates. 

First we discuss the simplest situation, where the 
ground state is non-degenerate, and later we shall gener- 
alize this to our situation with a degenerate ground state. 
The Chern number is in the context of quantum Hall sys- 
tems related to a measurable physical quantity, the Hall 
conductance. The boundary averaged Hall conductance 
for the (non-degenerate) ath many-body eigenstate of 
the Hamiltonian is [13, [1^: crf^ = C{a)e'^/h, where the 
Chern number C(a) is given by 



C{a) = r dO, r de2{d,At^ - d2A["^), (12) 

Jo Jo 



where -4j"' (6*1 ,6*2) is defined as a vector field based on the 
eigenstate ^''^"''(^i, 62) on the boundary torus Si x Si by 



^(")(01,02)=Z(VI/(")|A|VI/W 



(13) 



It should be noted that the wave function ^'(")(6'i,6i2) 
is defined up to a phase factor on the boundary-phase 
space. Therefore, ^'^"'(^'1, ^2) and e'-'^f^i'^^^^f^H^'i, ^2) 
are physically equivalent for any smooth function 

f {01,62). Under this phase change, ^j"^(6'i, 6*2) trans- 
forms like a gauge: 



Afi 



P2; 

h,62) 



Af {01,92) -d,f {61, 



(14) 



(11) 



Hence, the Chern number integral is conserved un- 
der this gauge transformation and it encapsulates general 
properties of the system. Chern numbers has been used 
extensively in the quantum Hall literature, for character- 
ization of the localized and extended states (Ref. [s^ and 
refs. therein). In this paper, we use the Chern number as 
an indicator of order in the system. Moreover, it enables 
us to characterize the ground state in different regimes, 
especially where the calculation of the overlap with the 
Laughlin wave function fails to give a conclusive answer. 
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FIG. 4: (a) Twist angles of the toroidal boundary condition. 
In addition to the flux going through the surface there may 
also be a flux inside the torus or going through the hole in the 
middle. When encircling these fluxes the wave function ac- 
quire an extra phase represented by the boundary conditions 
in Eq. (b) Redefining the vector potential around the sin- 
gularities: Aj is not well-defined everywhere on the torus of 
the boundary condition. Therefore, another vector field A' j 
with different definition should be introduced around each sin- 
gularity {01 ,62) of Aj. Aj and A' j are related to each other 
by a gauge transformation x a^nd the Chern number depends 
only on the loop integrals of x around those singularities re- 
gions, c.f., Eq. (fT5)) . 



Before explaining the method for calculating the Chern 
number, we clarify some issues related to the degeneracy 
of the ground state. In some systems, the ground state 
can be degenerate, this can be intrinsic or it can be as a 
result of the topology of the system. If the ground state is 
degenerate, we should generalize the simple integral form 
of Eq. to take into account the redundancy inside 
the ground state manifold. For example, in the case of 
a two-fold degenerate ground state, an extra gauge free- 
dom related to the relative phase between two ground 
states, and this freedom should be fixed. In other words, 
as we change the twist angles, we can not keep track of 
the evolution of both states, since one can not distinguish 
them from each other. Therefore, to uniquely determine 
the Chern number of the ground state(s), we should re- 
solve this gauge invariance, which is treated in Section 
O and [331 



It important to note that the degeneracy in the non- 
interacting regime is fundamentally different from the de- 
generacy in the interacting case. In the non-interacting 
limit, the degeneracy can be lifted by a local perturba- 
tion e.g. impurity, while in the hardcore case, the de- 
generacy remains in the thermodynamic limit 37]. The 
latter degeneracy in the ground state is a consequence 
of the global non-trivial properties of the manifold on 
which the particles move rather than related to a sym- 
metry breaking which happens in conventional models 
e.g. Ising model. The topological degeneracy is not a 
consequence of breaking of any symmetry only in the 
ground state, instead it is the entire Hilbert space which 
is split into disconnected pieces not related by any sym- 
metry transformation. With a general argument. Wen 
(38t showed that if the degeneracy of a chiral spin sys- 
tem moving on a simple torus is fc, then it should be fcf 
on a torus with g handles (Riemann surface with genus 
g), therefore the topological detergency is an intrinsic 
feature of the system. In particular, in the context of 
quantum Hall effect, this multicomponent feature of the 
ground state on a torus has a physical significance, while 
the single component ground state on a sphere bound- 
ary condition gives zero conductance, the torus geome- 
try with multicomponent ground state results in a correct 
conductance measured in the experiment, since the torus 
boundary condition is more relevant to the experiment. 
Changing twists angles of the boundary will rotate these 
components i nto each other and gives an overall non-zero 
conductance 



31 



As studied in a recent work by M. Oshikawa and T. 
Senthil [3^, as a universal feature, it has been shown 
that in presence of a gap, there is a direct connection 
between the fractionalization and the topological order. 
More precisely, once a system has some quasiparticles 
with fractional statistics, a topological degeneracy should 
occur, which indicates the presence of a topological or- 
der. Therefore, the amount of the degeneracy is related 
to the statistics of the fractionalized quasiparticles e.g. 
in the case oi v = 1/2, the two-fold degeneracy is related 
to 1/2 anyonic statistics of the corresponding quasiparti- 
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cles. Chern number has been also studies for spin 1/2 sys- 
tem on a honeycomb lattice 40] for identifying Abelian 
and non-Abelian anyons. Bellissard et al. [4l| studies 
Chern number for disordered Fermi system using a non- 
commutative geometry. 

To resolve the extra gauge freedom related to the two 
degenerate ground states, we consider two possibilities: 
I. lifting the degeneracy by adding some impurities, II. 
fixing the relative phase between the two states in the 
ground state. Below, we explore both cases. In the first 
case, we introduce some fixed impurities to lift the de- 
generacy in the ground state for all values of the twist 
angles. This is an artifact of the finite size of the system 
which we take advantage of. In the presence of pertur- 
bation, we show that the system has a topological order 
in spite of poor overlap with the Laughlin state. In the 
second approach, we use a scheme recently proposed by 
Hatsugai [i^, which is a generalized form for degen- 
erate manifolds. 



3.2. Resolving the degeneracy by adding impurities 

In this section, we introduce some perturbation in the 
finite system in form of local potentials (similar to local 
impurities in electronic systems) to split the degeneracy 
of the ground state and resolve the corresponding gauge 
invariance, which allows us to compute the Chern num- 
ber. Furthermore, the fact that we can still uniquely de- 
termine the Chern number in the presence of impurities, 
shows that the system retains its topological order, even 
when the impurities distort its ground state wavefunction 
away from the Laughlin wavefunction. 

In the context of the quantum Hall effect, the conven- 
tional numerical calculation of various physical quantities 
such as the energy spectrum, screening charge density 
profile, wave functions overlaps, and the density-density 
correlation functions, can not be used for understanding 
the transport properties of electrons in the presence of 
impurities (although useful for studying of isolated im- 
purities 0i[S|)- Recently, D. N. Sheng et al. ^ cal- 
culated the Chern number as an unambiguous way to 
distinguish insulating and current carrying states in the 
J/ = 1/3 quantum Hall regime which correspond to zero 
and non-zero Chern number, respectively. In this work, a 
weak random disorder was used to lift the degeneracy of 
the ground state (three- fold for i/ = 1 /3) for a finite num- 
ber of electrons. The energy splitting between the lowest 
three states then decreased with increasing number of 
particles, which indicates the recovery of the degeneracy 
in the thermodynamic limit. Moreover, the mobility gap 
can be determined by looking at the energy at which the 
Chern number drops towards zero. This energy gap is 
comparable to the energy gap obtained from the experi- 
ment and it is not necessarily equal to the spectrum gap 
which separates the degenerate ground state from the ex- 
cited states This shows the significance of Chern number 
calculations for understanding these systems. 



In a finite system, the coupling to a single-body inter- 
action, e.g. impurities, can lift the degeneracy and one 
can uniquely determine the Chern number for the indi- 
vidual states by direct integration of Eq. (H^). On one 
hand, the impurity should be strong enough to split the 
energy in the ground state (in this case E2,Ei, where Ej 
denotes the energy of the jth energy level) for all values 
of the twist angles. On the other hand, these impuri- 
ties should be weak enough so that the energy splitting 
in the ground state remains smaller than the thermody- 
namic gap (in this case E2 — Ei E^ — E^)- 

To calculate the Chern number of individual level, as 
mentioned in the pervious section, wc have to fix the 
phase of the wavefunction. The method that we explore 
in this section can be considered as a simplified version 
of the general method developed by Hatsugai ^] which 
we will explore in the next section. Following Kohmoto's 
procedure [i^, we assume that the ground state ^'(^i, 6*2) 
may be expanded for all twist angles on a s-dimensional 
Hilbert discrete space ^{9x,92) — (ci, C2, c^). If 
Aj{9i, O2) in Eq. is a periodic function on the torus 
of the boundary condition, then by application of Stoke's 
theorem, the Chern number will be always zeros. The 
non-triviality (non-zero conductance in the case of quan- 
tum Hall system) occurs because of the zeros of the wave 
function, where the phase is not well-defined. Therefore, 
-4(^1 J ^2) is not well defined everywhere and its flux inte- 
gral can be non-zero. To uniquely determine the Chern 
number, we assume that the wave function and the vec- 
tor field are not defined for certain points (6'", 62) in Sn 
regions on the torus of the boundary condition. For sim- 
plicity, we first discuss this procedure, in the case of a 
non-degenerate ground state. For calculating the inte- 
gral, we should acquire another gauge convention for the 
wave function inside these Sn regions, e.g. in a discrete 
system, we may require an arbitrary element of the wave 
function to be always real, and thereby we can define a 
new vector field -4'^"' {61 ,62), which is well defined inside 
these regions. These two vector fields differ from each 
other by a gauge transformation (Fig. 2]) : 

Afie,,92)-Af\e,,92)^d,xi9i,02), (15) 

and the Chern number reduces to the winding number 
of the gauge transformation x(^i;^2) over small loops 
encircling {0^,02), i.e. dSn, 

C(«) = E^/ ^X-de. (16) 

The one-dimensional gauge Eq. ([TS]) should be re- 
solved by making two conventions. For example, in one 
convention the first element and in the other the sec- 
ond element of the wavefunction in the Hilbert space 
should be real i.e. transforming the ground state ^ 
into *<j, = P$ = v]/\l/t$ where $ = (1, 0, O)'^ is 
a s-dimensional vector and P is a projection into the 
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ground state and similarly with the other reference vec- 
tor $' — (0, 1, 0)'''. Since the gauge that relates two 
vector fields is the same as the one that relates the cor- 
responding wavefunctions (similar to Eq. (jl5p ). we can 
uniquely determine the gauge transformation function x 
by evaluating ^{61, 82) = e'^ = $tp$'. Therefore, the 
Chern number will be equal to the number of windings 
of X around regions where = ^^p^ = j^-j^p jg ^eio. 
Counting the vorticities has a vigorous computational ad- 
vantage over the conventional method of direct integra- 



tion of Eq. (fT2|) . In the direct integration, we need to 
choose a large number of mesh points for the boundary 
angles, because of the discrete approximation of deriva- 
tives in Eq. , and this makes the calculation computa- 
tionally heavy. We note that for the system on a lattice, 
we should exactly diagonalize the Hamiltonian which is 
a sparse matrix as opposed to the continuum case where 
the Hamiltonian is a dense matrix residing on a smaller 
projected Hilbert space (lowest Landau level). 
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FIG. 5: (color online) Chern number associated to low-lying energy states, in the presences of impurities. Due to the impurity 
potential (a), two- fold degenerate ground state splits and the wavefunction overlap with the Laughlin state drops to 52% and 
65%, for the first and the second energy state, respectively. The results are for 3 atoms on a 6x6 lattice (a = 0.17) in the 
hard-core limit, (b) D,(Q\,Q-i) for the first level has no vorticity. However, for the second level, as shown in (c), n(Si,02) has 
vorticity equal to one associated with regions where either or A^ vanishes. 



For removing the degeneracy, in our numerical sim- 
ulations, we add a small amount of impurity which is 
modeled as delta function potentials located at random 
sites with a random strength, of the order of the tunnel- 
ing energy J. This is described by a Hamiltonian of the 
form R = UJ^iii., where i numerates the lattice site, 
hi is the atom number operator and f7™ is the strength 
of the impurity at site i. 

We choose reference states $ and $' to be eigenvectors 



of the numerically diagonalized Hamiltonian at two dif- 
ferent twist angles. In Fig. [51 vorticities of e*-^ associated 
with the first and the second energy level is depicted. It 
is easy to see that the Chern number associated to the 
two ground states is one. Number of vortices may vary 
for the first and second ground states, but their sum is 
always equal to one. The hard-core limit ([/ ^ J) is 
very similar to the case of fractional quantum Hall ef- 
fect, which in the context of Hall systems, means a share 
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of 1/2 (in /h unit) for each ground state [35|. When 
the onsite interaction strength is smaU {U < J), the ther- 
modynamic gap becomes comparable to the ground state 
energy sphtting E2 — Ei ^ — E2, the Chern number 
can not be uniquely determined, and the system doesn't 
have topological order. On the other hand, in the limit 
of strong interaction {U ^ J), the total Chern number 
associated to the ground states is equal to one, regardless 
of the impurity configuration. Moreover, in the hard-core 
limit, although the ground state is not described by the 
Laughlin wavefunction, since it is distorted by the impu- 
rity ( in our model it can be as low as 50%), the Chern 
number is unique and robust. This is an indication that 
the topological order is not related to the symmetries of 
the Hamiltonian and it is robust against arbitrary local 
perturbations 37]. These results indicate the existence 
of a topological order in the system and robustness of the 
ground states against local perturbations. 



3.3. Gauge fixing 

The method developed in the previous section has the 
graphical vortex representation for the Chern number 
which makes it computationally advantageous compared 
to the direct integration of Eq. (|12p . It can not, how- 
ever, be applied directly to a degenerate ground state, 
and therefore we had to introduce an impurity potential 
which lifted the degeneracy. On the other hand, a signif- 
icant amount of impurity in the system may distort the 
energy spectrum, so that the underlying physical proper- 
ties of the lattice and fluxes could be confounded by the 
artifacts due to the impurities, especially for large a. To 
address this issues, in this section, we explore a gener- 
alized method of the previous section based on Refs. li^ 
and[43l which works for a degenerate ground state. 

By generalizing the Chern number formalism for a de- 
generate ground state manifold, instead of having a sin- 
gle vector field ^j"^ (6*1, 6*2), a tensor field ^j"'''^ (6*1, 6*2) 



should be defined, where a,P — 1,2,. 
degenerate ground state 



, q for a q-fold 



A 



(«,/9) 



(^1 



d 



(vl;(")|_|vl;(« 



(17) 



Similar to the non-degenerate case, when W'"'^^ is not 
defined, a new gauge convention should be acquired for 
the regions with singularities. This gives rise to a tensor 
gauge transformation x'"''^''(^ii ^2) on the border of these 
regions 



A 



i)-A' 



^) = a,x("-«((?i,02). (18) 



Following Hatsugai's proposalf43| for fixing the ground 
state manifold gauge, we take two reference multiplets $ 
and $' which are two arbitrary s x q matrices; q is the 
ground state degeneracy (equal to 2 in our case). In our 
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FIG. 6; (color online) (a) shows the argument of 0(^1,612) 
as arrows for fixed $ and (b) and (c): surface plots of 
detA,j and detA^, (blue is lower than red). 61 and 62 changes 
from zero to 2n. These plots have been produced for 3 atoms 
with = 6 (a = 0.24) in the hard-core limit on a 5x5 lattice. 
The total vorticity corresponding to each of the reference wave 
functions ($ or "!>') indicates a Chern number equal to one. 
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FIG. 7: (color online) n(6'i,6'2) for fixed $ and <!>'. Oi and 
02 changes form zero to 2tt. This plot has been produced for 
4 atoms with Nij) = 8 in the hard-core limit on a 5x5 lattice 
{a — 0.32). Although, there are more vortices here compared 
to Fig. |6l the total vorticity corresponding to each of the trial 
functions ($ or $') indicates a Chern number equal to one. 

numerical simulation, we choose the multiplets to be two 
sets of ground state at two different twist angles far from 
each other, e.g. (0,0) and (tt, tt). We define an overlap 
matrix as = <!>''' P<I> where P = ^'^'^ is again the 
projection into the ground state multiplet, and consider 
the regions where detA$ or detA$/ vanishes (similar to 
zeros of the wave function in the non-degenerate case). 
Hence, the Chern number for q degenerate states, will be 
equal to the total winding number of Tr ;^("''^) for small 
neighborhoods, Sn, in which dctA$ vanishes 

C(l,2,...,(7) ^T.Y'f ^Tr x^"''') • (19) 

which is the same as the number of vortices of $') = 
det($'l'P$'). ft should be noted that the zeros of det A$ 
and detA'^ should not coincide in order to uniquely de- 
termine the total vorticity. In Fig.[6l we have plotted fl, 
detA$, and detA$/, found by numerical diagonalization 
of the Hamiltonian for a mesh (30 x 30) of winding an- 
gles 9i and 62- In this figure, the Chern number can be 
determined be counting the number of vortices and it is 
readily seen that the winding number is equal to one for 
the corresponding zeros of detA$ (or detA$'). 

We have calculated the Chern number for fixed 1^—1/2 
and different a's by the approach described above. The 
result is shown in Table HI For low a ^ 1 , we know from 
Sec. [2]that the ground state is the Laughlin state and we 
expect to get a Chern number equal to one. For higher 
a, the lattice structure becomes more apparent and the 
overlap with the Laughlin state decreases. However, in 
our calculation, the ground state remains two-fold degen- 
erate and it turns out that the ground state Chern num- 
ber tends to remain equal to one before reaching some 
critical ac — 0.4. Hence, also in this regime we expect to 
have similar topological order and fractional statistics of 



Atoms 


Lattice 


a 


Chern/state 


Overlap 


3 


6x6 


.17 


1/2 


0.99 


4 


6x6 


.22 


1/2 


0.98 


3 


5x5 


.24 


1/2 


0.98 


3 


4x5 


.3 


1/2 


0.91 


4 


5x5 


.32 


1/2 


0.78 


3 


4x4 


.375 


1/2 


0.29 



TABLE I: Chern Number for different configurations in the 
hard-core limit for fixed filling factor u = 1/2. The Laughlin 
state overlap is shown in the last column, although it deviates 
from the Laughlin state. Although the ground state deviates 
from the Laughlin state, the Chern number remains equal to 
one half per state before reaching some critical Qc — 0.4 where 
the energy gap vanishes. 



the excitations above these states on the lattice. 

For the arguments above to be applicable, it is essen- 
tial that we can uniquely identify a two-fold degenerate 
ground state which is well separated from all higher ly- 
ing states. For higher flux densities, a > ac, the two-fold 
ground state degeneracy is no longer valid everywhere 
on the torus of the boundary condition. In this regime, 
the issue of degeneracy is more subtle, and the finite size 
effect becomes significant. The translational symmetry 
argument [2^ , which was used in the Section [21 is not 
applicable on a lattice and as pointed out by Kol et al. 
[47| the degeneracy of the ground state may vary period- 
ically with the system size. Some of the gaps which ap- 
pear in the calculation may be due to the finite size and 
vanish in the thermodynamic limit, whereas others may 
represent real energy gaps which are still present in the 
thermodynamic limit. To investigate this, we study the 
ground state degeneracy as a function of boundary an- 
gles (01, 6*2) which are not physical observable and there- 
fore the degeneracy in thermodynamic limit should not 
depend on their value. In particular. Fig. [5] shows the 
energy levels of five particles at a = 0.4 for different val- 
ues of the twist angles. The first and the second level 
are split at {9i = 6*2 = 0), while they touch each other 
at (6*1 =02 = tt). We have observed similar behavior for 
different number of particles and lattice sizes e.g. 3 and 
4 atoms at a = 0.5. In this case, the system seems to not 
have a two-fold degeneracy. Therefore, the ground state 
enters a different regime which is a subject for further 
investigation. 

For having the topological order, it is not necessary 
to be in the hard-core limit. Even at finite interaction 
strength U ^ Ja, we have observed the same topological 
order with the help of the Chern number calculation. If 
U gets further smaller, the energy gap above the ground 
state diminishes (as seen in Sec. and the topological 
order disappears. 

We conclude that the Chern number can be unambigu- 
ously calculated for the ground state of the system in a 
regime where Laughlin's description is not appropriate 
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FIG. 8; (color online) Low-lying energy levels as a function of 
twist angles. For high a the degeneracy of the ground state is 
a function of twist angles. The shown results are for 5 atoms 
on a 5x5 lattice i.e. a = 0.4 (a) shows first three energy 
manifolds as a function of the toroidal boundary condition 
angles, (b) shows a cross section of (a) at 62 — n for seven 
lowest energy levels. The first and the second energy levels 
get close to each other at 9\ = 62 = tt. 



for the lattice. The non-zero Chern number of a two-fold 
degenerate ground state, in this case equal to one half 
per state, is a direct indication of the topological order 
of the system. 



4. EXTENSION OF THE MODEL 

In Sections [2] and [3] above, we have investigated the 
conditions under which the fractional quantum Hall ef- 
fect may be realized for particles on a lattice. The mo- 
tivation for this study is the possibility to generate the 
quantum Hall effect with ultra cold atoms in an opti- 
cal lattice but the results of these sections are applicable 
regardless of the method used to attain this situation. 



In this and the following sections, we investigate some 
questions which are of particular relevance to ultra cold 
atoms in an optical lattice. First, we introduce a long 
range, e.g., dipole-dipole, interaction which turn out to 
increase the energy gap and thereby stabilizes the quan- 
tum Hall states. We then turn to the case of v = 1/4 
and show that in order to realize this state, it is essential 
to have some kind of long range interaction. 



4.1. Effect of the long-range interaction 

In an experimental realization of the quantum Hall ef- 
fect on a lattice, it is desirable to have as large an en- 
ergy gap as possible in order to be insensitive to exter- 
nal perturbations. So far, we have studied effect of the 
short range interaction, and we have shown that the gap 
increases with increasing interaction strength, but the 
value of the gap saturates when the short range interac- 
tion becomes comparable to the tunneling energy J . 

In this section, we explore the possibility of increasing 
the gap by adding a long-range repulsive dipole-dipole 
interaction to the system. Previously, such dipole-dipole 
interaction has also been studied in Ref. 
onlinecitecooper2005 as a method to achieve Read-Rezayi 
states [4^ of rapidly rotating cold trapped atoms for 

= 3/2 and as a means to realize fractional Quan- 
tum Hall physics with Fermi gases [i^. The dipole- 
dipole (magnetic or electric) interaction is described by 
the Hamiltonian: 



d-d 



Udipoh 



E 

l<i<j<N 



p, • pj - 3(ny • p,)(nij • Pj) 



where n.y = (r^ - rj)/|ri - Vj 



(20) 

Pi are unit vectors repre- 



senting the permanent dipole moments and the position 
vectors are in units of the lattice spacing a. For sim- 
plicity, we assume that all dipoles are polarized in the 
direction perpendicular to the plane. With time indepen- 
dent dipoles, the strength of the interaction is given by 



U dipole 



(or 



r) where /i's (p's) are the perma- 



nent magnetic (electric) dipole moment. Static dipoles 
will thus give repulsive interaction Udipoie > 0, but ex- 
perimentally time varying fields may be introduced which 
effectively change the sign of the interaction [sO]. For 
completeness, we shall therefore both investigate posi- 
tive and negative Udipoie, but the repulsive interaction 
corresponding to static dipoles will be the most desirable 
situation since it stabilizes the quantum Hall states. 

Experimentally the dipole-dipole interaction will nat- 
urally be present in the recently realized Bose-Einstein 
condensation of Chromium [sij which has a large mag- 
netic moment. However, for a lattice realization, polar 
molecules which have strong permanent electric dipole 
moments is a more promising candidate. For typical po- 
lar molecules with the electric moment p 1 Debye, on 
a lattice with spacing a 0.5/zm, Udipoie can be up to 
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FIG. 9; (color online) (a) The overlap of the ground state 
with the Laughlin wavefunction (dashed lines) and four low 
lying energies of the system (solid lines) versus the dipole- 
dipole interaction for four atoms on a 6x6 lattice, (b) Gap 
enhancement for a fixed repulsive dipole-dipole interaction 
strength Udipoie ~ 5J versus a. The results are shown for 
iV = 2(n), N = 3(*), iV = 4(+) 



few kHz, an order of magnitude greater than the typical 
tunneling J/inh which can be few hundreds of Hz J^J. 

To study the effect of the dipole-dipole interaction, we 
again numerically diagonalizc the Hamiltonian for a few 
hardcore bosons [U ^ J), in the dilute regime a ^ 0.3, 
while varying strength of the dipole-dipole interaction. 
The results of this simulation is shown in Fig. [5] (a). 
After adding the dipole interaction, the ground state in 
a dilute lattice remains two-fold degenerate, since the in- 
teraction only depends on the relative distance of the par- 
ticles and keeps the center of mass motion intact. If the 
dipole interaction becomes too strong {\Jdipoie 3> J) , the 
ground state wave function deviates from the Laughlin 
wavefunction, but the topological order remains the same 
as for the system without dipole interaction. We veri- 
fied this by calculating the Chern number as explained 
in Sec. [3] for different values of the dipole-dipole interac- 
tion strength \J dipole and observed that the total Chern 
number of the two-fold degenerate ground state is equal 
to one. Moreover, as it is shown in Fig. [H] (b) adding 
such an interaction can increase the gap: the lower curve 
corresponds to the hard-core limit discussed in the pre- 
vious work (To| and the upper curve corresponds to the 
system including the dipole-dipole interaction. This en- 



hancement varies linearly with the flux density a in a 
dilute lattice and doesn't depend on the number of par- 
ticles and consequently, it is expected to behave similarly 
in the thermodynamic limit. 

One of the impediment of the experimental realization 
of Quantum Hall state is the smallness of the gap which 
can be improved by adding dipole-dipole interaction. In 
this section, we showed that this improvement is possible 
and moreover, by Chern number evaluation, we verified 
that adding dipole interaction doesn't change the topo- 
logical behavior of the ground state manifold. 



4.2. The case oi v = 1/4 

So far we have concentrated on the case oiv— 1/2. In 
this section, we briefly investigate the case oi v = 1/4. 
It is expected that the Laughlin wavefunction remains a 
good description for the ground state of a bosonic sys- 
tem for any even g, where v = 1/q. Following Haldane's 
argument [2^, due to the center of mass motion, the de- 
generacy of the ground state is expected to be g-fold on 
a torus. Similar to the case of = 1/2, the Laughlin 
wavefunction should be a suitable description for any q 
provided that the magnetic field is weak so that we are 
close to the continuum limit, i.e. a ^ 1. Also the Chern 
number is expected to be equal to one for the q-fold de- 
generate ground state, which in the context of quantum 
Hall effect means a share of 1/g of the conduction quan- 
tum e'^/h for each state in the q-fold degenerate ground 
state manifold. 

We have done both overlap and the Chern number cal- 
culation to check those premises. In the case of = 1/4, 
significant overlap occurs at low a < 0.1. The average 
wave function overlap of four lowest energy eigenstates 
with the Laughlin wavefunction is depicted in figure [TUl 
where we have used a generalization of the Laughlin 
wavefunction for periodic boundary conditions similar to 
Eq. O [13. 

We observe that the Laughlin wavefunction is a reliable 
description of the system with = 1/4 but only for much 
more dilute lattices {a < 0.1) compared to = 1/2 where 
significant overlap occurs for a < 0.3. Contrary to v = 
1/2, where the gap is a fraction of the tunneling energy 
J, the gap for i/ ~ 1/4 between the 4- fold degenerate 
ground state and the next energy level is infinitesimal. 
The reason for the vanishing gap can be understood in 
the continuum limit from the argument put forward in 
Ref. [5^ : as noted previously the Laughlin wavefunction 
is an exact eigenstate of the Hamiltonian with an energy 
per particle equal to the lowest Landau level energy. The 
energy of the in — A state is thus equal to the m — 2 state. 
It thus costs a negligible energy to compress the = 1/4 
state to the v — 1/2 state, and therefore there is no gap. 
In an external trap the system will always choose the 
highest density state which is the v — \/2 state. Note, 
however that this argument only applies to short range 
interactions. For long range interactions, we expect to 
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FIG. 10: (color online) The overlap of the first four low- 
lying energy states with the Laughlin wavefunction for the 
case of = 1/4 on a torus. The dashed (dotted) line shows 
the overlap for the system without (with) dipole interaction 
{Udipoie = 5J). The Laughlin state is only a good descrip- 
tion for a more dilute lattice a < 0.1 compared to v = 1/2. 
The dipole interaction stabilizes the system and the overlap 
is more significant for higher values of a < 0.2. 



see a non vanishing gap. 

Even though that with short range interactions, the 
gap is very small in our numerical calculations, it is still 
sufficiently large that it allows us to unambiguously de- 
termine the Chern number for the ground state mani- 
fold as described in Sec. 13.31 As expected the calculation 
shows that the Chern number is equal to one correspond- 
ing to a four-fold degenerate ground state consistent with 
the generalization of the fcrmionic case in the fractional 
quantum Hall theory [31, |35[. In Fig. [TOl the overlap 
of the first four lowest energy state with the Laughlin 
wavefunction is depicted. In the absence of the dipole in- 
teraction, the ground state overlap is significant only for 
a ^ 0.1, however, by adding a moderate dipole interac- 
tion {Udipoie — 5J), the overlap becomes more significant 
for a larger interval of the flux density, i.e. a < 0.25. This 
is due to the fact that state with lower density become 
more favorable in the presence of a long-range repulsive 
interaction. 

We observed that adding a dipole interaction would 
lead to an improvement of the gap for v = 1/2 and make 
the Laughlin overlap more significant for larger interval of 
magnetic field strength a in the case oi i' — 1/4. There- 
fore this long-range interaction can be used as an tool for 
stabilizing the ground state and make the realization of 
these quantum states, experimentally more feasible. 



and imaging the momentum distribution. In Ref. 1^, it 
was shown that this technique provide some indication 
of the dynamics in the system. This measurement tech- 
nique, however, only provides limited information, since 
it only measures the single particle density matrix, and 
provides no information about the correlations between 
the particles. In Refs. [s^ and[5^ a more advanced mea- 
surement techniques were proposed, where the particle 
correlation is obtained by looking at the correlations in 
the expansion images. In this section, we study Bragg 
scattering as an alternative measurement strategy which 
reveals the excitation spectrum of the quantum system. 
In Ref. [H bosonic quantum Hall system responses to a 
perturbative potential is studied. We focus on Bragg 
scattering where two momentum states of the same in- 
ternal ground state are connected by a stimulated two- 
photon process [5^] . By setting up two laser beams with 
frequencies cui and uj2 and wave vectors fci and k2 in the 
plane of the original lattice, a running optical superlat- 
tice will be produced with frequency loi — uj2 and wave 
vector ki — k2- (Both frequencies lji and lj2 should be 
close to an internal electronic dipole transition in the 
atoms). The beams should be weak and sufficiently de- 
tuned so that direct photon transitions are negligible, i.e. 
El, £2, 7 ^ — Wo, UJ2 — where is the frequency 
of the transition, 7 is the corresponding spontaneous de- 
cay rate and £1 , £2 are the Rabi frequencies related to the 
laser-atom coupling. In this perturbative regime, the in- 
elastic scattering of photons will be suppressed; the atom 
starts from the internal ground state, absorbs one pho- 
ton from e.g. beam 1 by going to a virtual excited state 
and then emits another photon into beam 2 by return- 
ing to its internal ground state. After this process, the 
system has acquired an energy equal to hw = hiuji — UJ2) 
and a momentum kick equal to q = ki — k2- There- 
fore, the overall effect of the recoil process is a moving 
AC Stark shift as a perturbing moving potential, and 
the effective Hamiltonian represents the exchange of the 
two-photon recoil momentum and the energy difference 
to the system and is proportional to the local density i.e. 
H cx p(r)e~'('^i-'^2)*+*('="i-'=''2)-^-h c.c. 

This process can be used to probe density fluctuations 
of the system and thus to measure directly the dynamic 
structure factor S{q,uj) and static structure factor S{q). 
This kind of spectroscopy has been studied for a BEC in 
a trap by Blakie et al. [SQ] and Zambelli et al. l57il and 
has been realized experimentally in Refs. [H, l59l. l60l. l6ll . 



5. DETECTION OF THE QUANTUM HALL 
STATE 

In an experimental realization of the quantum Hall 
states, it is essential to have an experimental probe which 
can verify that the desired states were produced. In most 
experiments with cold trapped atoms, the state of the 
system is probed by releasing the atoms from the trap 



Also, Bragg spectroscopy in an optical lattice is dis- 
cussed in Ref. |6J in a mean-field approach and also in 
Ref. 65i as a probe of the Mott-insulator excitation spec- 
trum. On the other hand, in the context of quantum 
Hall effect, the static structure factor has been studied 
for probing the magnetoroton excitations [6^ and charge 
density waves [67| . The dynamic and static structure 
factors are given respectively as 
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FIG. 11: (color online) (a) Structure factor and (b) energy 
spectrum for a 11x11 lattice with 3 atoms on a torus. Points 
shows the momentums allowed by the boundary conditions. 
The dotted line in (b) shows for comparison the low energy 
spectrum of a free particle which equals J{qaf' . 



lattice sites. Fig. [TlTa) shows the structure factor for 
the case v — 1/2 ior a. small a calculated from our nu- 
merical diagonalizations. In the data presented here, we 
have chosen Ij = qx and but the result should be simi- 
lar in other directions in the lattice plane. We see that 
S{q) is modulated at a momentum corresponding to the 
magnetic length. For the parameters that we have in- 
vestigated, the general features of the structure factor is 
independent of the size of the system. 

We obtain the excitation spectrum shown in Fig. [11] 
(b) similar to Ref. [6^ by the Feynman-Bijl approxima- 
tion. In the continuum limit {a <ti 1), we assume that 
the low-lying excitations are long-wavelength density os- 
cillations and their wave functions can be approximated 
to have the form cx pfc|0). Therefore, the variational es- 
timate for the excitation energy is u}{q) ~ fi^g^/ 2mS{q). 
At zero momentum and at the momentum corresponding 
to the magnetic length, there are gaps, and we also ob- 
serve a deviation from the free particle spectrum similar 
to the magneoroton case as a reminiscent of the Wigner 
crystal. It should be noted that the deviation does not 
depend on the size or the number of particles in the sys- 
tem. As clearly seen in the Fig. [TT]the energy spectrum 
and structure factor deviate from those of free particles, 
therefore, it could be used as an experimental probe of 
the system. 

The structure factor and excitation spectrum imply 
some general features that are very different from that 
of the Mott-insulator and superfluid states, and can be 
used a powerful experimental indication of the quantum 
Hall states. 



GENERATING MAGNETIC HAMILTONIAN 
FOR NEUTRAL ATOMS ON A LATTICE 



S{q,uj) 



E 

n,0 



\{n\p\q)\Q)\H{uj~E,, + Eo), (21) 



s{q) = Y.\^n\pHm? = E K^^i E ^''-^mw m) 

n,0 n,0 fi 

where the density fluctuation operator is defined as 
pH'I) ^ X^m n -^<j(™j ")4nCn and the coefficient are de- 
fined as Fourier transforms of the Wannier functions: 
Aq{m,n) — J (Pr e^^'^(l)*{r ~ rjn)4>{r — fn), where the 
Wannier function (/'(r — r„) is the wave function of an 
atom localized on a site centered at r„. Below, we focus 
on deep optical lattices, where Aq{m,n) = e*'' '"™^™.™- 

In the structure factor, there is a sum over the ex- 
cited states \n) and ground states |0) and the self-term 
is thus excluded. The ground state on a torus is two- 
fold degenerate and therefore in our numerics, we add 
the contribution of both. 

Since we are working on a discrete lattice, there will 
be a cut-off in the allowed momentum given by the lat- 
tice spacing (/max = 7r/a, where a is the distance between 



Recently, there has been several proposals for produc- 
ing artificial magnetic field for neutral atoms in opti- 
cal lattices [l3, H^, however, the implementation of 
each of them is still experimentally demanding. Recently, 
there has been an experimental demonstration of a rotat- 
ing optical lattice |23j which is equivalent to an effective 
magnetic field (see below). The technique used in this 
experiment, however, generates a lattice with a large lat- 
tice spacing, because it uses laser beam which are not ex- 
actly counter propagating. This longer spacing reduces 
the energy scale in the lattice and thereby also reduces 
quantities such as energy gaps. Here, we shall now in- 
troduce an alternative method for generating a rotating 
optical lattice, which does not result in an increased lat- 
tice spacing. This method consists of rotating the optical 
lattice by manipulating laser beams. 

In a frame of reference rotating with angular velocity 
uj around the z — axis, the Hamiltonian for a particle of 
mass m in an (planar) harmonic trap of natural frequency 
ujo is 
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— 7r/4 < {0 = Lot) < 7r/4, where the parameters are de- 
fined in Fig. [T^ If we consider a square lattice with 



FIG. 12: Proposal for realizing a rotating optical lattice. Four 
AOMs (black boxes) changes the direction of the lattice beam, 
which are subsequently focussed in the middle of the setup by 
four lenses (grey). Simultaneously varying the four diffraction 
angles in the AOMs will generate a rotating optical lattice. 



H = 



(p — muz X r)^ 
2m 



ujz.r X p 



V 
1 

—'i 
2 



At resonance loq = lo the form is equivalent to the Hamil- 
tonian of a particle of charge q experiencing an effective 
magnetic field B = sj x (mojzxr/q) — {2muj/q)z. There- 
fore, by simply rotating the optical lattice, we can mimic 
the magnetic field for neutral atoms. 

To rotate the optical lattice, we propose to set up four 
acousto-optic modulators (AOM) and four focusing com- 
posite lenses as shown in the Fig. [121 By sweeping the 
acoustic wave frequency, the beams can be focused and 
make a rotating optical lattice. 

In an AOM, for the first order diffracted light we have 
sin 9b = ^ , where A is the wavelength of sound in 
the medium, A is the light wavelength and 6b is half 
of the angle between a diffracted beam and the non- 
diffracted beam (Fig. [121) . By increasing the frequency of 
the acoustic wave, the diffraction angle increases. How- 
ever, the beam should be focused by a large aperture lens 
so that it always passes the region where we want to make 
the optical lattice. By focusing a similar but counter- 
propagating beam, we can make a rotating standing wave 
(Fig. [T2I) . By repeating the same configuration in the 
transverse direction, we can make a rotating optical lat- 
tice. In particular, if the AOM is far from the composite 
lenses, D ^ d, then x/D = A/ A and x/d = tan 9 where 



dimensions Nx 



given by rotation is 



TV, the number of magnetic flux 



BA _T:Af^ 
$0 ~ 2 uj, 



(24) 



where ujr = hk"^ /2M is the atomic recoil frequency. 

On the other hand, the upper limit for the magnetic 
field comes from the single Bloch band approximation 
which we made in writing the Hamiltonian for the opti- 
cal lattice. In order for the particles to remain in the first 
band, the traveling lattice beams should move them adi- 
abatically. From Ref. l68l, the adiabaticity condition for a 
moving lattice with an acceleration 77 equal to uP'MX/^ 

at the edge, is mryA <C ft^, where LOp is the frequency 
difference between the first and the second band in the 
lattice. This puts a limit on how large the lattice can 

become N <^ — . 

Hence, for = 1/2 with a lattice filling fraction ^ 
I and a typical recoil frequency ujr — (27r)4 kHz, one 
can enter the regime of fractional quantum Hall effect by 
rotating the lattice at w ~ (27r)650 Hz. If a deep optical 
lattice is used e.g. ujp ~ lOwr, the adiabaticity condition 
is easily satisfied for a lattice of size N ^ 1000. The 
experimentally challenging part will, however, likely be 
to mitigate the instability of the lattice caused by the 
thickness of the beam and aberration of the lenses at the 
turning points i.e. near 9 — 7r/4. 



7. CONCLUSIONS 

An extended study of the realization of the fractional 
quantum Hall states in optical lattices has been pre- 
sented. We showed that a Hamiltonian similar to that of 
a magnetic field for charged particles can be constructed 
for neutral ultra cold atoms and molecules confined in 
an optical lattice. By adding an onsite interaction for 
the case oi v — 1/2, an energy gap develops between 
the ground state and the first excited state, which in- 
creases linearly as all and saturates to its value in the 
hardcore limit U ^ J. We learned that the Laughlin 
wavefunction is a reliable description of the system for 
low flux densities a < 0.25. However, for higher a's, the 
lattice structure becomes more pronounced and a bet- 
ter description of the system can be carried out by in- 
vestigating the Chern number associated to the ground 
state manifold. The Chern number indicates that the sys- 
tem has topological order up to some critical flux density 
ac — 0.4, where the properties of the ground state man- 
ifold starts to change. We have also studied v — 1/4, 
where compared to v = 1/2, the Laughlin wavefunction 
only describes the ground state for lower values of the flux 
a < 0.1. We showed that a dipole-dipole interaction can 
enhance the gap and stabilize the system, and therefore 
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make the ground state more experimentally realizable. 
Bragg spectroscopy has been studied as a potential ex- 
perimental tool to diagnose theses incompressible states. 

Characterization of the ground state by evaluating the 
Chern number, developed in Sec. 13. 3i can be general- 
ized to other interesting many-body systems where the 
conventional overlap calculation fails to work. In partic- 
ular, this method can be applied to ground states with 
non-Abelian statistics which are appealing candidates for 
fault-tolerant quantum computation. 
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